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We describe the conditional and unconditional dynamics of two coupled quantum dots when 
one dot is subjected to a measurement of its occupation number using a single electron transistor 
(SET). The measurement is made when the bare tunneling rate though the SET is changed by the 
occupation number of one of the dots. We show that there is a difference between the time scale for 
the measurement-induced decoherence between the localized states of the dots and the time scale 
on which the system becomes localized due to the measurement. A comparison between theory and 
current experiments is made. 



I. INTRODUCTION 

There have recently been a number of suggestions for a quantum computer architecture that use quantum dots of 
varying kinds If these schemes are to be practical many important physical questions need to be answered, one 

of which is how to readout physical properties such as charge or spin at a single electron level QJj| . In this paper we 
present a quantum trajectory analysis of a general scheme to readout a single electronic qubit using a single electron 
transistor (SET). We adopt a general phenomenological description of the SET in which the tunneling rate through 
the SET is conditioned on the occupation or otherwise of a nearby quantum dot. 

We consider two spatially separated quantum dots which are strongly coupled so that delocalized states of their 
relevant degrees of freedom can form. To be specific, we imagine each dot to have a single electronic bound state 
that can be occupied. Thus the average occupation number of each dot must be less than unity. This restriction can 
easily be removed to account for spin, or multiple electron states. We label each dot with an index 1, 2 and let Cj, cj 
represent the Fermi annihilation and creation operators for each single electron state (see figure 1). 
The two dots are strongly coupled via the tunnel coupling Hamilton 

V = ih^(c{c 2 - 4d) (1) 
Thus the total Hamiltonian of the two-dot system is 

2 

H = h^Ui4ci + V (2) 

i=l 

In what follows we will work in an interaction picture and assume that the energies of each bound state are equal 
(again this can be relaxed). Coulomb blockade effects have been ignored at this stage, but can easily be included 
without significantly changing the results of this paper. 

The single particle eigenstates of this Hamiltonian are even and odd superpositions of the bare states of each well. 
Such states are thus delocalized over the two-dot system and are sometimes called 'molecular states' in the literature 
||. The localized states can then be represented as an even and odd superposition of the delocalized states. The 
localized states are not stationary; rather, the system will periodically oscillate between them. That is to say the 
system will tunnel coherently between the two dots. 

To this coherent system we add a measurement device which determines the presence of an electron on one of the 
dots, say dot 1, which we shall refer to as the target (see figure [l]). The model is based on a SET tunnel junction 
containing a single bound state on the island. The interaction between the target and the SET is via a Coulomb 
blockade. Thus the interaction Hamiltonian between the SET and the target must commute with the target electron 
number operator , c\c\. This makes it a QND (quantum nondemolition measurement) of electron number |Q. The 
Coulomb blockade changes the current flowing through the SET. In simple terms if there is no electron on the target 
the island state is biased so as to allow little or no current to flow though the SET. This is the quiescent state of the 
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SET. However when there is an electron on the target, the coulomb blockade shifts the bound state on the island to 
allow a greater current to flow through the device (see figure ||) 

We derive a master equation to describe the behaviour of the target system. This master equation describes the 
unconditional evolution of the measured system when the results of all measurement records (that is current records) 
are averaged over. This will tell us the rate at which coherence in the target system is destroyed by the measurement. 
However we also need to know how the system state depends on the actual current through the device in order to 
determine how quickly the conditional state of the electron becomes localized which is measure of the quality of the 
measurement. 

One approach to this problem is to keep track of many different states of the system, corresponding to the different 
numbers of electrons which have tunneled through the SET. This is the approach used for example in Ref . || . Here 
we adopt an alternate method which gives a more intuitive picture for the conditional dynamics. We use a conditional 
stochastic master equation which gives the evolution of the measured system, conditioned on a particular realization 
of the measured current. The instantaneous state of the target conditions the measured current while the measured 
current itself conditions the future evolution of the measured system in a self consistent manner. This approach to 
measurements has been variously called the quantum trajectory method |^] or quantum monte carlo method |lo| . 



II. SET MODEL 

Consider a two dot system with the coupling in Eq. ([!]). If the electron is in dot-2 the quiescent rate of current 
tunneling though the SET is a constant which we will denote Dq. However if there is an electron on dot-1, the rate 
of tunneling through the SET changes to Do + D\ with D\ > 0. If Dq docs not equal zero then a current spike 
(resulting from a tunnel event in the SET) does not necessarily imply that the electron in the measured system is in 
dot-1. In an ideal device the quiescent tunneling rate, Dq is zero. In reality Johnson noise on the circuit containing 
the SET will give a non-zero quiescent tunneling current. 

Assuming that the SET island state can be adiabatically eliminated, it is possible to derive a master equation for 
the state of the coupled dot system. This is done in the appendix, and the result is 

% = -i[ H > P\ + ldccT>[c\c 1 }p = Cp (3) 
where the irreversible part is defined for arbitrary operators A and B by 

V[A]B = J[A]B - A[A]B, (4) 



where 



The decoherence rate is given by 



J\A\B = ABA"* (5) 
A[A]B = -(A^AB + BA i A). (6) 



lAec = 2D +D 1 . (7) 



The fact that the irreversible term is a function of the number operator in the target qubit is an indication that this 
describes a QND measurement of the occupation of the dot-1. It is easy to verify that the stationary solution of this 
master equation is an equal mixture of the two accessible electronic states. 

The stochastic record of measurement ideally comprises a sequence of times, being the times at which electrons 
tunneled through the SET. In practice of course these events are not seen due to a finite frequency response of the 
circuit (including the SET) which averages each event over some time. However for the purpose of this paper we will 
take the zero response time limit. In this limit the current consists of a sequence of 5 function spikes. Formally we 
can write i(t) = edN/dt, where dN(t) is a classical point process which represents the number (either zero or one) 
of tunneling events seen in an infinitesimal time dt, and e is the electronic charge. We can think of dN(t) as the 
increment in the number of electrons N(t) in the collector in time dt. It is this variable, the accumulated electron 
number transmitted by the SET, which is used in Ref. ||. 

The point process dN(t) is formally defined by the conditions 

[dN{t)} 2 = dN(t) (8) 
E[dN(t)]/dt = A)Tr[(l - m)p c (i)(l - m)] + (D + Di)Tt[mp c (t)ni] 

= D + D 1 (m) c (t). (9) 
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Here E[x] denotes a classical average of a classical stochastic process x, and 

ni = cjci (10) 

is the occupation number operator for the first dot. The first of these equations simply expresses the fact that dN(t) 
equals zero or one. The second says that the rate of events is equal to a background rate D plus an additional rate 
D\ if and only if the electron is in the first dot. 

In Eq. (||), the system state matrix p c (t) is not the solution of the master equation (|^). That is because if one has 
a record of the current dN/dt through the SET then one knows more about the system than the master equation 
indicates. That is to say, p c (t) is actually conditioned by dN(t') for t' < t, hence the subscript c. The first way of 
writing Eq. (|^) hints at how dN(t') conditions p c (t). From the appendix, the state at time t + dt given dN(t) = 1 is 

p x {t + dt) = dt [ d (i - m)p(t)(i - m) + (d + D 1 )n l p{t)m] (11) 

This is an unnormalized state whose norm is equal to the probability of that event (dN{t) = 1) occurring, as seen 
above in Eq. (^|). 

The normalized state can be written more elegantly as 

, , , {Do + DiJ[ni] + 2D V[n 1 ])p(t) _ (D J[1 - m] + (Dp + D^Jjn^pjt) 

Tr{(D + D 1 J[n 1 } + 2D V[ni])p(t)} E[dN(t)/dt] ' 1 ' 

where J is as defined above in Eq. 

To write the full conditioned evolution we need to know the density operator po(t + dt) given that dN(t) — 0. This 
can be found from Eq. (|l^) plus the fact that when averaged over the observed classical point process dN, 

p (t + dt) + pi{t + dt) = (1 + Cdt)p(t). (13) 

That is to say, on average the system still obeys the master equation (||) . From this equation, we obtain 

p (t + dt) = P (t) - dt{D A[l - m]p(t) + {D + D 1 )A[n 1 )p{t) + i[H, p(t)}}, (14) 
= p(t) - dt {D oP (t) + D 1 A[n 1 }p{t) - i[H, p(t)]} , (15) 

where A is as defined above in Eq. (^J) . Once again this state is unnormalized and its norm gives the probability that 
dN(t) = 0, that is 

Tr[p (t + dt)} = 1 - E[dN(t)] (16) 
Using the variable dN(t) explicitly, the conditioned state at time t + dt is 

^^--^M^ +[i - im] M^k- (17) 

Since dN(t) is almost always zero we can set dN(t)dt — and expand this expression to finally obtain the stochastic 
master equation, conditioned on the observed event in time dt 



dp c = dN(t) 



Do + D 1 J[n 1 } + 2D V[n 1 } _ ' 
D + DiTr[p c ni] 



Pc 



+ dt + DiTrip.n^p ~ i[H, p c }} (18) 

Note that averaging this equation over the observed stochastic process (by setting dN equal to its expected value) 
gives the unconditional master equation (0). 



III. AVERAGE STEADY STATE PROPERTIES 



We now analyze in some detail the ensemble averaged properties of the system based on the unconditional master 
equation. In particular we calculate the stationary noise power spectrum of the current through the SET when there 
is the possibility of coherent tunneling between dot-2 and the measured dot-1. The details of how the quantum 
stochastic processes in the SET determine the average current though the SET are given in reference |L2| . The link 
with the stochastic formalism of the preceding section is that the current i(t) through the SET is given by 
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(19) 



First we calculate the steady state current 

ioo = E[£(t)]oo 

= e(D + D 1 (c\c 1 } oo ) (20) 
= e(D + ^-), (21) 

where the 00 subscript indicates that the system is at steady-state. The fluctuations in the observed current, i(t) are 
quantified by the two-time correlation function: 

G(r) = E[i(t)i(t + T)-i 2 00 } oc (22) 
= ei oo 8(r) + (i(t),i(t + T))£ . (23) 

Here (A, B) = (AB) — (A) (B). The fact that the multiplier of the shot noise is eioo rather than the usual (e/2)ioo 
is because of the approximation we have made in treating the SET. Specifically, we have adiabatically eliminated the 
SET by taking the limit where any electron which tunnels onto the SET island from the emitter immediately tunnels 
off to the collector. This means that, on the time scales we are interested in, there is a perfect correlation between 
the emitter current and collector current. This leads to a doubling of the shot noise level. Of course, at very high 
frequencies, higher than we are interested in, the true shot noise level of (e/2)zoo could still be seen in principle. 

To relate these classical averages to the fundamental quantum processes occurring in the well we apply the theory of 
open quantum system PJTl[| to the present system. Specifically, we can relate the correlation function for the current 
to the following quantum averages 

(i(t)i(t + r)}^° = e 2 Tr [(£>„ + D x J[m] + 2D V[n 1 ])e CT (D + D x J\n x \ + 2D V[n 1 ])p oo ] . (24) 

Because Poc is an equal mixture of the two electron states, it satisfies 

V[ ni } Poc = 0. (25) 

In addition, the following identies for arbitrary operators A and B are easy to prove: Tr [P[/l]£>] = 0, Tr [j7[ni]_B] = 
Tr[mB], Tr[e CT B] = Tr[B], and Tr[Ae CT Poo ] = Tr[A Poo }. Using these simplfications we obtain 

(i(t),i(t + T)) T *° = Die 2 {Tr[nie CT J[ni} Poo \ -Ti[n lPoo } 2 } . (26) 
Evaluating this expression we find 

G{t) = eiocdir) + — — ===- , (27) 

8 \ V(7dec/4) 2 - O 2 J 

where 

P± = -( 7d ec/4) ± ^{ ldcc /±y-n\ (28) 

and where the first term represents the shot noise component as discussed above. The power spectrum of the noise is 

S{u)= drG(r)2cos(wr), (29) 
Jo 

which evaluates to 



e 2 Dln 2 /2 
V(7dec/4) 2 - fi 2 P 2 



= zioo + „ )J = \ -a— a - -3—2 \ ■ (30) 



In the case that f2 > 7dec/4 the spectrum will have a double peak structure indicating that coherent tunneling is 
taking place between the two coupled dots. For smaller only a single peak appears in the spectrum. We can thus 
use the noise power spectrum of the current though the SET as a means to measure the tunnel coupling between dots 
if the tunnel coupling is high enough. We illustrate this in figure [| 
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IV. ANALYTICAL RESULTS FOR CONDITIONAL DYNAMICS 



We now return to the stochastic master equation for the conditioned state, 

' D + DiJicld] +2D V[c\ Cl } 



dp c = dN 



Pc 

1 (31) 



D Q + DxTr[p c c\ci] 
+ dt |--Di^{ctci,p } + DiTripcclc^p - i[H,p c ]^ . 

Comparing this to the unconditional master equation 

p = -i[H, p] + ^ dcc V[c\ci]p (32) 

we see that decoherence between the two coupled dots, 1 and 2, takes place at the rate 7d oc = 2D + £>i, but that 
the system decides between the two possibilities (electron on dot-1 or on dot-2) on a time scale that depends on D\ 
and Dq in some more complicated way. Of course this measurement time scale is necessarily at least as large as 
the decoherence time scale because successfully distinguishing between the two dots would by definition destroy any 
coherence between them. 

The different measurement time scales can be derived most easily by introducing the Bloch representation of the 
state matrix: 

p = - (I + xa x + yay + za z ) (33) 
where the Pauli matrices are defined using the Fermi operators for the two dots 

<?x = cjc 2 + c\ci (34) 
<jy = —ic\c2 + ic\c\ (35) 
a z = c\c 2 - c\c x (36) 

In this representation the means of the Pauli matrices a a are given by the respective coefficient a, with a = x, y, z. 

The stochastic master equation can now be written as a set of coupled stochastic differential equations for the Bloch 
sphere variables as 

dZc = n Xc dt + ^(i - zl)dt - mt) D ^-^ z ] )l2 07) 

dx c = —ilz c dt ^-z c x c dt — dN(t)x c (38) 

dy c = ~^-z c y c dt ~ dN(t)y c . (39) 

Again the c-subscript is to emphasize that these variables refer to the conditional state. If we average over the noise, 
the ensemble dynamics is then seen to be given by 

dz „ , , 

- = ftc (40) 

dX Tdcc / j -i \ 

*=- nz -— x (41) 

± = -l^y (42) 
dt 2 y v ' 

where a = E[a c ] denotes the averaging over the ensemble of conditional states. These equations are exactly what 
would be obtained directly from the ensemble averaged master equation Eq(^|) . In particular we note that the average 
population difference z between the dots is a constant of the motion in the absence of any free Hamiltonian. However 
the stochastic differential equations enable us to calculate important averages that are not obtainable from the master 
equation. For example, if the model does indeed describe a measurement of c\c\ = (1 — cr z )/2, then, in the absence of 
tunneling, we would expect to see the conditional state become localized at either z = 1 or z = — 1. Indeed for O = 
we can see from the conditional equation for z c that z c = ±1 is a fixed point. 
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We can take into account both fixed points by considering z 2 . In the absence of tunneling this must must approach 
1 for all trajectories, since the system will eventually become localized due to the measurement in one dot or the 
other. Therefore it is sensible to take the ensemble average E[z 2 ] and find the rate at which this deterministic quantity 
approaches one. Noting that for a stochastic variable we have d(z 2 ) = 2zdz + dzdz, and that E[d/V 2 ] = E[cL/V] = 
[D Q + Di(l - z c )/2}dt, we find that 



dt 



= E 



D 2 (l-z 2 ) 2 



4£>o + 2 J D 1 (l 



(43) 



If the system starts state which has an equal probability for single electron to be on each dot then z c (0) = and in 
the ensemble average this would remain the case. However if we ensemble average z 1 over many quantum trajectories 
then for short times we find 

E[z 2 (5t)} = St (44) 

That is to say, the system tends towards a definite state (with z c = ±1 so z 2 = 1) at an initial rate of Z? 2 /(4Z?o + 2£>i). 
For vanishing Dq, this is the same as the decoherence rate, D\j2, as expected. But for Dq 3> D\, the rate goes to 
(Di/2Dq) x Di/2 <C Dx/2. That is, the rate at which the system becomes localized at one or the other dot is much 
less than the decoherence rate. This result cannot be obtained from the ensemble averaged master equation alone. 
It is a direct reflection of the fact that for Dq / a tunneling event cannot be unambiguously attributed to the 
location of an electron on the double dot system. As the rate of localization is a direct indication of the quality of 
the measurement, we can use the localization rate defined as 

7loc = 4^T2A (45) 

as an important parameter defining the quality of the measurement. This parameter is related to the signal-to-noise 
ratio for this measurement as we now show. 

For a Poisson process at rate R, the probability for to events to occur in time T is 

p( m; T) = ffi^ e -* T (46) 
to! 

The mean and variance of this distribution are equal and given by E(m) = Var(TO) = RT . Now consider an electron 
which is, with equal likelihood, in either dot, so that z = 0. If the electron is in dot 1 then the rate of electrons 
passing through the SET is Dq + D\\ if it is in dot two then it is just Dq. These two possibilities will begin to be 
distinguishable when the difference in the means of the two distributions pirn, T) is of order the square root of the 
sum of the variances. That is, when 



D\T ~ sj D Q T + (D + Di)T (47) 
Solving this for T gives a characteristic rate 

D 2 

T- 1 =-i . (48) 

The right hand side of this expression is simply twice the 7i oc defined above. A similar conclusion is reached in 
reference ||. 

In the ideal limit of no quiescent current in the SET Do = 0, the stochastic master equation can be replaced by 
a stochastic Schrodinger equation, and will collapse to a single possibility at a rate D\ which is the same as the 
decoherence rate. The effect of Dq is most clearly seen in the other limit, D\ <C Dq, as noted above. In this limit 
the single electron makes only a small relative change in the tunneling rate through the SET. As the rate of jumps 
also becomes large, the trajectories in this limit take on the appearance of diffusion rather than jumps. The rate at 
which the electron localizes into one well or the other scales as D~^ 2 Dq, which is much longer than the decoherence 
time scale -D^ 1 . 
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V. NUMERICAL SIMULATIONS OF CONDITIONAL DYNAMICS 



We now turn to numerical simulations of the conditional evolution and to estimate the conditions for a good 
measurement. Unlike traditional condensed matter measurements we wish to describe repeated measurements made 
on a single quantum system rather than a single measurement made upon an ensemble of systems. To do this we use 
the conditional dynamics of the measured system given a particular measurement record as described by the above 
stochastic evolution equation ( |3l| ) 

We return to the Bloch description defined in Eq. (33). In what follows we will assume that y(0) = 0. For the form 
of tunneling used here the value of y does not in fact change under either conditional or ensemble averaged dynamics. 
If the conditional state of the system remains in a pure state then x\ + yl + z% — x\ + z 2 = 1 . As noted previously 
this can only occur if the bare tunneling rate (-Do) is zero, when a tunneling event can unambiguously be attributed 
to the occupation of the target dot and no information is lost about the state of the system. In the more realistic case 
in which Dq ^ 0, we can use the quantity x\ + z c 2 as a measure of the purity of the sate, or equivalently as a measure 
of how much information the conditional record of measurements gives about the actual state of the two coupled dots. 
If the conditional state is a maximally mixed state of a two state system then x\ + z 2 = 0. We now describe in detail 
the numerical simulation of the conditional dynamics. 



A. No Background Current 

First we consider the case Dq — 0, so the system is always in a pure state. Typical trajectories are shown in 
Fig. 4 for various values of fl. For small <C -Di/2 we see little evidence for coherent tunneling. Most of the 
time the electron is localized almost entirely in one well or the other. However, there is an asymmetry between the 
wells. A transition from dot-2 into dot-1 (the target) is sudden, occurring whenever an electron tunnels through the 
SET. A transition the other way takes a time of order 2/ 'D\. This time is still much smaller than the average time 
between state-changing transitions, which can be shown analytically to be Di/Q 2 . Thus over a long time, as shown 
in Fig. 4(a), the system still has the appearance of a random telegraph. This behaviour gives rise to a single-peaked 
noise spectrum, as shown in Fig. 3(a). 

For moderate fl ~ D\/2 the system is no longer well- localized in one dot or the other. Rather, the dynamics is 
complicated with clearly non-sinusoidal oscillations from one dot to the other interspersed with jumps into the target 
dot. The fact that oscillations are present gives peaks in the current noise spectrum, as shown in Fig. 3(b). The 
position of these peaks is a frequency less than fi, as shown analytically in Sec. III. 

For f2 3> Di/2 the dynamics once again becomes simple, with nearly sinusoidal oscillations interspersed with jumps 
which occur with an average rate of D\j2. This corresponds to a noise spectrum having a very sharp feature at 
u> m ±fi, as shown in Fig. 3(c). 

The change in behaviour as fl increases is summarized in Fig. 5. There we plot E[z^] versus Q/D\. The quantity 
E[z^] measures how well localized the electron is at one well or the other, and would be 1 if the electron were always 
localized and if it were never localized. It is actually possible to calculate this quantity numerically without using 
a stochastic ensemble, as follows. 

With no background current, every time an electron tunnels through the SET the electron on the dots is known to 
be on dot 1. If there are no further SET tunneling events for a time t later then from Eq. (|l4|), the system evolves up 
to that time by the equation 

dp (t) = -dt{D 1 {4c Ll po(t)}/2 + i[H,p a (t)]}. (49) 

Because there is no background current, and because the initial state is pure, it is possible to rewrite this in terms of 
a non-Hermitian Schrodinger equation 

d\Mt)) = -dt(iH + D lC \ Cl /2)\Mt)). (50) 

Here it must be remembered that the norm of this state represents the probability for the event that no electron has 
passed through the SET since the last one a time t ago: 

Po(t) = Tr[p (t)] = (Mt)\Mt))- (51) 



It is not difficult to show that the solution to Eq. (50) satisfying the initial condition V>o(0)) = 1 1 , 0) is 



\Mt)) =a(t)|l,0)+/?(t)|0,l), (52) 
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where the occupation numbers refer to the dots one and two in order. Here a and (3 are real numbers defined by 

1 



a(t) 



A+ - A- 



( A+e A + f - A_e A - 4 ) 



-fi/2 



A + - A_ 



where 



A+ = 




(53) 
(54) 

(55) 



The conditioned quantum expectation value for a z is 

(Mt)\a z \Mt)) _ 1 



z (t) 



Po(t) 



(56) 



Now in steady state the probability po(t) that there is no increment in N(t) for a time t ago is related to the 
probability qo(t) that the last increment was a time t ago by 



qo (t) 



Po(t) 



J °°p (s)ds 



(57) 



Since at steady state all conditioned states are uniquely identified by how long it has been since the last SET event, 
the ensemble average for z 2 is simply given by 



EUtf = / q (t)[z (t)} 2 dt 
IQ 



(f3 2 + a 2 ) dt 



f3 2 +a 2 



-dt. 



(58) 
(59) 



Unfortunately it does not appear possible to evaluate the second integral here analytically. However a numerical 
integration is easy. The results, shown in Fig. 5, is in agreement with the ensemble averages obtained numerically 
using the stochastic master equation. 



B. A Finite Background Current 



We next consider the case where Dq ^ 0. We show two plots, both with Q = D%, which is a regime in which 
coherent tunneling is clearly evident in the current noise spectrum. The first plot, in Fig. 6, is for Do = D\. Here 
coherent oscillations are still evident in z, but z rarely attains its extreme values of ±1. The conditioned state is no 
longer pure, even immediately after a count. Also, the conditioned state following a count now depends on the state 
before the count. For this reason an exact solution by the method of the preceding section is impossible. 

The second plot, in Fig. 7, is for Dq = lO-Dj.. Here coherent oscillations are no longer obvious in the condition 
mean of cr 2 , even thought they are present in the spectrum (as small features above the shot noise), as calculated in 
Sec. III. In this regime Dq 3> D\ so the diffusive limit discussed at the end of Sec. IV applies. 



VI. DISCUSSION AND CONCLUSION 



The three parameters we need to compare our theoretical results with experiment are f2, Dq and D\. The two 
incoherent tunneling rates can be obtained by considering how they determine the steady state current through the 
device. This is given in Eq (|2l|). Our model implicitly assumes that the quiescent noise in the SET is shot-noise 
limited, based as it is on elementary tunneling events. However from the point of view of the macroscopic circuit in 
which the SET is placed, the tunnel junctions appear as a capacitor in series with a resistor ( see note 14 in reference 
J13|). If the resistance of the junction is R and the capacitance is C the fundamental time constant for the junction is 
■j^j. This sets an upper bound for the tunneling rate 7 < -^j For a typical Al/AlO^/Al junctions we have C « 0.2/F 
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and R ~ 50fcil and thus the time constant is 7 w 10 11 s -1 . For a double tunnel junction device, the maximum 
conductance is achieved for a symmetric pair of tunnel junctions. The value of D\ in this case is given by D± = 7/2 
where 7 is the tunneling rate of the SET under conditions of maximum conductance (see appendix). We thus estimate 
that Di « 5 X 10 10 s _1 . The background tunneling rate through the SET depends on temperature as well as the bias 
conditions (see appendix). Typical maximum and minimum conductance for the SET at different temperatures have 
been measured by Joyez et al. Jig] . At a temperature of lOOmK the minimum conductance is approximately zero, 
and thus at this temperature we can safely take Do = 0, the zero temperature result. However at a temperature of 
400mK the ratio of the maximum to minimum conductance is 2.2. This indicates that at 400 mK, Do = 4x 10 10 s _1 . 

The value we choose for the tunneling rate depends strongly on the particular quantum dot system. We will consider 
the value appropriate for the single electron measurement scheme of Kane et al. |Q . In this model the two localized 
states correspond to an electron on a tellurium ion donor or at a nearby interface below the donor. The tunneling 
rate is then between the donor state and the interface state. Kane et al. estimate that for this system f2 w 10 9 s _1 . 
This value is too low to observe a peak in the noise power spectrum away from zero frequency. However the tunneling 
rate could be increased by changing the donor interface bias voltage. 

In conclusion we have presented a simple model to describe single electron measurements on a coupled double 
quantum dot system using an SET. We have given both the stationary (ensemble averaged) properties of the current 
through the SET as well as the conditional dynamics of repeated measurements on a single system. This illustrates 
how quantum trajectory methods may be naturally adapted to single electronics, and aid in the interpretation of 
ensemble averaged properties. We believe these models will prove useful in current attempts to fabricate quantum 
logic gates in solid state devices. 



APPENDIX A: DERIVATION OF THE MASTER EQUATION. 

It will suffice to consider a single quantum dot near the SET. This allows us to remove any reference to the electron 
field labeled by Ca, c\. The SET is modeled as a single biased double barrier (single well) device with a single bound 
state on the well described by the Fermi operators b, b< . The total Hamiltonian for the system including the reservoirs 
is 

H = H + H CB + H RT + H LT (Al) 

The term Hcb is the Coulomb blockade term and is given by || 

Hcb = h X c{c x tfb (A2) 

where fix IS the Coulomb blockade energy gap (see figure |^) . Note this term can only be nonzero if there is an electron 
on the island and on the dot, in which case the energy of the island electron is increased. The terms H R t, Hlt 
described the tunneling between the many modes in the left and right ohmic contacts and the bound state on the 
SET H 

k 
k 

where a^ k , a Rk are respectively the Fermi field annihilation operators for the left and right reservoir states at 
momentum k. The tunneling matrix elements between respectively the left and right Ohmic contacts and the island 
are T^ k , T Rk . The free Hamiltonian for the the system is 

H = h*^2 w fe°ife a £Ai + u k a Rk a Rk + Tiujic\ci + Tiuj tfb (A3) 

We now transform to an interaction picture to remove the terms Hq + Hcb- The dynamics in the Schrodinger picture 
is now described by the time dependent Hamiltonian 

k 

+T* Lk a Lk tfe- lxtc i Cl e l ^- Ua]t 

T Rk a Rk be lxt ^ Cl e- l{ ^^ a)t 

+T R _ k a Rk tfe- txtc i Cl e i{uJ "^ o)t 
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Using the fact that (cjci)™ = c\c\ we find 

H I {t)=H 1 {t) + H 2 {t) 

where 

H x {t) = (1 - c\ Cl ) Y.iTLka^be-' 1 ^-^ + H.c.) + (T Rk a Rk be' l ^^)t + jj. c .) 
fc 

iJ 2 (i) = cjd ^(Tifco^fte-*^-" -^* + H.c.) + (T^a^fee" 1 ^-^-^* + H.c.) 
fe 

Notice that if there is no electron on the dot and c\ci — > then the second term is zero and the first term is a 
standard tunneling interaction onto a bound state with energy fitOQ. On the other hand if there is an electron on the 
dot c\c\ — * and the first term is zero and the second term is a standard tunneling interaction onto a bound state 
with energy Ti{ujq + \) as expected. 

The derivation of the master equation for the state matrix R for the system (SET and quantum dot) can now 
proceed using standard techniques which we will sketch. The objective is to obtain a semigroup evolution in Lindblad 
form (that is to say positivity-preserving irreversible dynamics) for the state of the SET island and the dot alone 
with no reference to the ohmic contacts. The ohmic contacts are treated as perfect Fermi thermal reservoirs with 
a very fast relaxation constants. Each ohmic contact (left and right) remains in thermal equilibrium with chemical 
potentials pl, PR, but the total system is not in thermal equilibrium due to the external bias potential, V with 
eV = pi — PR (see references Jl^JlJ f° r further discussion). We first define a time interval 5t which is slow compared 
to the dynamics of the island and the dot but very long compared to the time scale in which the ohmic contacts relax 
back to their steady state. The change in the state matrix W of the system (SET and dot) and environment (Ohmic 
contacts) from time t to t + St, to second order in the tunnel coupling energy, is given by 

rt + 6t 

W{t + 6t) = W{t)-i6t[H I (t),W(t)]-6t dti[fr/(i),[£fi(ti),W(ti)]] (A4) 



We now make the first Markov approximation and assume that at any time the state of the total system may be 
approximated by W(t) = R(t) ® pl® Pr , that is to say the left and right ohmic contacts instantaneously relax back 
to Fermi distributions. We now obtain an evolution equation for R(t), the state of the island and the dot by tracing 
over the reservoirs. The result is 



dR(t) 
dt 



+ llLfM + Th/hK)] 2^(1 - c\ Cl )]R 

+ [7^(1 - f L (uo + X )) + 7a(1 - /h(wo + X))] V[bc\ Cl ]R 

+ [y L f L (Lu a + X ) + 7fl./fl(^o + X)} V[tfc\ Cl ]R 

where for arbitrary operators A and B, V[A]B = ABA^ — ^(A^AB — BA^A) and where /l,_r(w) is the Fermi filling- 
probability for the left/right ohmic contact at the energy tku. The rates Jl,r and "f' L R determine the rate of injection 
from the left ohmic contact into the island or emission from the island into the right ohmic contact under the conditions 
of no electron on the dot (unprimcd) and with an electron on the dot (primed) . These are evaluated using the second 
markov approximation as 

lL = \T Lko \\ (A5) 
lR = \T Rko \ 2 , (A6) 
j'L = \T Lk ' \ 2 (A7) 
i'r = \Tr K \ 2 (A8) 



where fco = y / 2mujo/h and k' = ^/2m(u>o + xj/ft. 

The ideal quiescent state of the SET is defined as /l(wo) — 1> /h( w o) = 1 while /z,(wo + x) = 1, /i?(^o + x) = 0- 
Under these conditions the master equation reduces to 

^ = lR V[b{\ -c\ Cl )]R + 7L X>[&t(l - c \ Cl )]R + i R V\bc\ Cl ]R + i L V[tfc\ Cl ]R (A9) 
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We now wish to derive a master equation for the state matrix p for the dot alone. This is easiest if we assume that 
IRiIr are much larger than all other rates in the system. In this case it is possible to adiabatically eliminate the 
SET island using techniques similar to that in Ref. [fyl. We expand the state matrix R in powers of I/7.R or l/~f' R as 

ie = po®|0}(0|+pi®|0)(0|. (A10) 

The equations of motion for pi and po are 

pi = -KrA[1 - n{\pi + jlJ[1 - m]po ~ j' R A[ni]pi + j' L J[ni]p (All) 
Po = 1rJ[1 - "i]Pi - j L A[l - rii]po + l R J[ni]pi - 7i^[ n i]Pi (A12) 

Here ri\ = c\c\ and J and A are as defined in Eqs. (Q), (^|). Under the above conditions, we can slave p\ to po so that 

{ 7r A[i - ni ] + 7R^[m])pi = (ilJ[i - m] + tlJMpo (A13) 

Operating on both sides alternately by ^7"[ni] and J7[l — tii] it is easy to show that 

i R J[nx]pi = i L J[n{\ P o (A14) 
7rJ[1 - ni]pi = ilJ[1 - ni]po. (A15) 

Substituting these into Eq. ( A12| ) yields 

p a = {i L +lL)V[ ni ] Po (A16) 

Since the probability or their being an electron on the SET is very small we can say p ~ po. Hence we have derived 
The master equation (||) (without the Hamiltonian term) for the dot alone 

p={2D + D 1 )V[c\c 1 ]p. (A17) 

Here we have defined 

0i=7t-7£, (A18) 

A) = 7l- (A19) 

Because the SET collector reservoir in the two cases (an electron on the dot and an electron not on the dot) are 
independent (due to the SET energy shift), the state conditioned on an electron entering the collector is an incoherent 
mixture of the two possible paths. From quantum trajectory theory P,|ll[, the unnormalized state conditioned on 
this event is 

dt (7 R J[b(l - ni)] + i R J[bni}) R. (A20) 

The norm of this state matrix gives the probability for this event, and is equal to the norm of 

dt{ lR J[l-n l ]+i R J[n 1 ]) Pl . (A21) 

^From the adiabatic elimination procedure above, this is equal to 

dt{7 L J[l-n 1 }+ 1 ' L J[n x ])p. (A22) 

This is the unnormalized state p\(t + dt) of the dot alone conditioned on an electron tunneling through the SET. 
From this it is easy to derive the rate of such tunnelings as 

7l (1 - ni> + i L (ni) - D + D x (c^ . (A23) 
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FIG. 1. Schematic representation of a single electron measurement for two coupled quantum dots 

FIG. 2. SET using coulomb blockade for a single electron measurement. The Coulomb blockade gap is labeled Ecb and 
the tunneling rates in the 'on' position are 7l and through the left and right barriers respectively. In (a) the electron is 
localized on dot-2 and a background current eDo flows. In (b) the electron is localized on dot-1 and the current e(Do + Di) 
flows in the SET. 

FIG. 3. A plot of the noise power spectrum normalized by the shot noise level for Do = 0. (a) f2 = 0.1, (b) f2 = 0.5, (c) 
Q = 5.0. 

FIG. 4. A plot of the conditional population difference dynamics of z c (t) versus scaled time for various values of the tunneling 
rate. In all cases Do = and time is measured in units of D^ 1 . (a) Q. — 0.1, (b) Q = 0.5, (c) fl = 5.0. 

FIG. 5. A plot of the average of the conditional quantity zj? versus fl/Di, wich measures the extent to which the measurement 
localises the state of the dot. We set Do = 0. The solid lines refer to the exact result Eq(59) 

FIG. 6. The effect of finite background current with Do = 1-0 and Q. = 1.0. A plot of (a) the 'purity measure' x\ + yl + zl 
versus scaled time (units of D]~ ) and (b) the conditional population difference z c (t) for a typical trajectory. 

FIG. 7. The effect of finite background current with Do = 10.0 and = 1.0. A plot of (a) the 'purity measure' x\ + y* + z\ 
versus scaled time (units of D^ 1 ) and (b) the conditional population difference z c (t) for a typical trajectory. 
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